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The effective string picture of confinement is used to derive theoretical predictions for the interquark potential 
at finite temperature. At short distances, the leading string correction to the Unear confining potential between 
a heavy quark-antiquark pair is the "Liischer term". We assume a Nambu-Goto effective string action, and work 
out subleading contributions in an analytical way. We discuss the contribution given by a possible "boundary 
term" in the effective action, comparing these predictions with results from simulations of lattice Z2 gauge theory 
in three dimensions, obtained with an algorithm which exploits the duality of the Z2 gauge model with the Ising 
spin model. 



1. EFFECTIVE STRING PICTURE 

The effective string picture is expected to pro- 
vide a good physical description of confining 
gauge theories in the low energy regime. Assum- 
ing that two confined color charges are joined by 
a thin flux tube that fluctuates like a vibrating 
string, and describing the string world sheet dy- 
namics by an effective action S'^^j, one can derive 
quantitative predictions about the potential be- 
tween a confined quark-antiquark pair. In par- 
ticular, for massless string fluctuations, the sim- 
plest choice is the Nambu-Goto string action: 
i^off = a ■ A, which is proportional to the area 
A of the world-sheet surface and a is the string 
tension, appearing as a parameter of the effective 
theory. 

We consider a 3D system with the extension 
X L with Lg >> L. With periodic boundary 
conditions employed (at least) in the short direc- 
tion, the temperature T is proportional to 1/L. 

Taking into account leading quantum fluctua- 
tions, the result for the expectation value of the 
Polyakov loop correlator is given by 
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where 77 is Dedekind's function. The term associ- 
ated with the minimal world sheet surface induces 
the exponential area-law fall off, and the conse- 
quent linear rise in the interquark potential V{R). 
The first non-trivial contribution in S^a results in 
the determinant of the Laplace operator, and the 
corresponding contribution to the interquark po- 
tential V{R) — in a regime of distances shorter 
is the Liischer term P^: 
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Further terms in the expansion of the Nambu- 
Goto action give rise to a contribution involving 
a combination of Eisenstein functions 3 : 
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However, it is worth mentioning that, while Eq. 
Q and Eq. ^ are generally accepted, on the 
other hand the further contribution to V{R) ex- 
pressed by Eq. ^ is still under debate. 

A possible "boundary term" in the effective ac- 
tion, as proposed in ref. can be treated by 
means of a perturbative expansion in h (a parame- 
ter proportional to the coefficient of such a term) , 
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which induces a leading order correction hke 0]: 
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with a short distance contribution to V{R): 
— 
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2. Z2 LATTICE GAUGE THEORY 

We run numerical simulations of Z2 lattice 
gauge theory in three space-time dimensions 0]. 
We chose this model because the effective string 
picture is believed to be independent of the un- 
derlying gauge group, and Z2 gauge group is in- 
teresting from the point of view of the center rele- 
vance in confinement. Moreover, the duality with 
the Ising spin model and a novel algorithm allow 
us to reach high precision results. 

The pure 3D lattice gauge model is described 
in terms of a^.f^ £ Z2 variables defined on the 
lattice bonds. The dynamics is governed by the 
standard Wilson action, which enjoys Z2 gauge 
invariance. The partition function is expressed in 
terms of plaquette variables and reads: 



Z(/3)=^e-'^^ = ^exp 
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There are different phases: a confined, strong cou- 
pling phase, with massive string fluctuations for 
j3 < 0.47542(f) 0, a confined, rough phase, with 
massless string fluctuations (this is the regime 
we studied in our simulations) and a deconfined 
phase for /3 > 0.76f4f34(2) (6). 

This model is dual with respect to the Z2 Ising 
spin model in 3D. The ratio of Polyakov-loop 
correlators G{R) in the gauge model can be ex- 
pressed as a ratio of partition functions Zl^r in 
the spin model: 
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The index L x R means that links perpendicular 
to an L X i? plane have an antiferromagnetic cou- 
pling constant. The ratio of partition functions 
can be expressed as an expectation value in one 



of the two ensembles. However the correspond- 
ing observable shows extremely large fluctuations. 
This problem can be overcome by a factorization 
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where -^[lx^rj+i+i h^-s just one more antiferro- 
magnetic link than .^[lxaJ+i- The expectation 
values corresponding to these ratios can be eas- 
ily obtained. A similar method (snake algorithm) 
has been employed in refs. |7] to compute the 't 
Hooft loop in the SU (2) gauge model. 

Important features of our implementation are 
multi-level updating and a hierarchical organiza- 
tion of sublattices. The CPU time is roughly 
proportional to the inverse temperature L, and 
independent of the distance R between the quark 
sources, thus the algorithm is particularly useful 
for large interquark distances. 

3. NUMERICAL RESULTS 

Let F{R, L) be the free energy of a heavy 
quark-antiquark pair at finite temperature: 
G{R) — e^^^^-'^\ We studied "quantum terms" 
in free energy differences, by measuring the quan- 
tity defined as: 



Q{R, L) = F{R + l,L)- F{R, L) - aL 



(8) 



Fig. ^shows that for large interquark distances 
(namely: L < 2R) our numerical results are in 
good agreement with the NLO prediction eq. © 
from Nambu-Goto string. It is important to note 
that the agreement between numerical results and 
the NLO prediction is not the result of a fitting 
procedure. A "purely classical" area law is defi- 
nitely ruled out, and the LO term alone is not suf- 
ficient to describe the data. We also found that 
the coefficient for a possible "boundary term" in 
the effective action for this model is very small, 
suggesting that it is exactly zero, i.e. 6 = 0. 

Fig. 12 — see 0] for details — shows the de- 
viation of a quantity related to the free energy 
second derivative from the free string behaviour 
in the L > 2R regime. The normalization is cho- 
sen so as to allow a meaningful comparison among 
different LGT's. Our Z2 results and SU{2) data 
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Figure 1. Q{R,L) for i = 10 (i.e. T = ATc/b) at 
13 = 0.75180. z is defined as: z = Solid curves 
correspond to the free bosonic string prediction, 
while dashed lines are the NLO Nambu-Goto cor- 
rection. A pure area-law would be described by 
Q = 0. 



are in good agreement: this may be a possible 
signature of the center relevance to confinement. 

4. CONCLUSIONS 

We tested the effective string predictions with 
precise numerical data for the finite tempera- 
ture Zi2 lattice gauge theory, using an algorithm 
that exploits the duality properties of the model. 
We explored a wide range of distances and de- 
tected next-to-leading order effects. Our numer- 
ical results seem to rule out a possible "bound- 
ary term" in the effective string action describ- 
ing the present model. For large distances R and 
high temperatures (i.e. small L) we see an ex- 
cellent agreement of the numerical data with the 
NLO prediction Eq. Q. In the regime of short 
distance and low temperature, we compared our 
results with different gauge theories, which have 
been studied by other authors. We focused onto 
the deviation of a well-suited combination of free 
energy differences with respect to the free string 
prediction. In particular, the agreement between 
our Z2 results and SU{2) data might be inter- 
preted as a possible signature of the fact that the 



Figure 2. Behaviour of SU{i) (2j (white squares), 
SU{2) (crosses), and Z2 0] (black squares) 
gauge models at short distances. 



center degrees of freedom play an important role 
in the confinement mechanism. 
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